Semiclassical theory of energy diffusive escape 
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Thermal escape out of a metastable well is considered in the weak friction regime, where the bot- 
tleneck for decay is energy diffusion, and at lower temperatures, where quantum tunneling becomes 
relevant. Within a systematic semiclassical formalism an extension of the classical diffusion equation 
is derived starting from a quantum mechanical master equation. In contrast to previous approaches 
finite barrier transmission also affects transition probabilities. The decay rate is obtained from 
the stationary non-equilibrium solution and captures the intimate interplay between thermal and 
quantum fluctuations above the crossover to the deep quantum regime. 



I. INTRODUCTION 



Thermally activated escape over high energy barriers 
can be found in a huge variety of physical and chemical 
processes [H, 0, [3, 3 ■ While these issues have been stud- 
ied for quite a while now, the subject has gained new 
interest in the context of quantum information process- 
ing, where read-out devices for superconducting quan- 
tum bits based on Josephson junctions (JJs) have been 
implemented [1, IH these realizations the expo- 

nential sensitivity of the switching rate out of the zero 
voltage state on the barrier topology has been exploited 
to discriminate the qubit states. Recently, a new type 
of detector (Josephson Bifurcation Amplifier) has been 
developed where an underdamped microwave driven JJ 
is operated in a bistable regime close to a bifurcation of 
two stable dynamical states [1, Q . In this range the J J 
acts as an anharmonic oscillator which is extremely sen- 
sitive to thermal and quantum fluctuations and may thus 
switch from one dynamical state (small-amplitude oscil- 
lation) to the other one (large-amplitude oscillation) or 
vice versa. Even though in this type of switching process 
the escape happens to occur across a dynamical rather 
than an energy barrier, in a moving frame picture rotat- 
ing with the frequency of the external driving force, the 
situation can be mapped on a time independent escape 
problem with a non-standard Hamiltonian [lo| . 

Accordingly, a classical rate theory has been devel- 
oped already two decades ago [ll| ■ It is founded on 
Kramers' seminal approach [12l | to calculate the decay 
rate from the steady solution of a corresponding time 
evolution equation for the probability density. This ap- 
proach has been substantially extended in the 1980s Q], 
in particular, to lower temperatures where below a so- 
called crossover temperature quantum tunneling prevails 
against thermal activation [ij, 13 ■ In between there is 
a broad range where both quantum and thermal fluctu- 
ations are strongly intermingled. The details of this be- 
havior depend very much on the friction strength though 
"^[IBl- Qualitatively, two regimes must be distinguished 
3|, namely, the range of moderate to strong friction, 
where classically the bottleneck for escape is spatial dif- 
fusion, and the range of weak friction, where the escape 
is controlled by energy diffusion. In the former one, 



typically realized in condensed phase systems, quantum 
corrections have been derived within dynamical formu- 
lations as well as within thermodynamical approaches 

More involved is the situation in the energy diffusive 
regime, found in gas phase systems like e.g. ensembles of 
atoms or molecules, since at lower temperatures an ap- 
propriate time evolution equation in energy space is not 
at hand a priori. In some approaches the classical energy 
diffusion equation has been extended ad hoc by adding a 
loss term which captures the finite transmission through 
the potential barrier near its top [l^, The contin- 
uum version of a quantum mechanical master equation 
has been the starting point for other studies [111, [23, [l^ , 
where transition matrix elements have been treated clas- 
sically, however, while finite transmission and reflection 
probabilities appear as quantum mechanical input. It 
was shown that quantum fluctuations lead to an intricate 
behavior for low friction and that in certain ranges of pa- 
rameter space the reflection near the barrier top may even 
prevail so that the escape rate shrinks below its classical 
value. Classical transition probabilities have also been 
used in [2^ and inserted in a discrete version of a master 
equation. In [lo| a master equation strictly valid only 
for harmonic systems has been applied to study quan- 
tum effects in the switching between dynamical states 
as described above. Hence, what is really missing, is a 
consistent semiclassical derivation of an energy diffusion 
equation for escape out of a metastable well in the regime 
of intermediate temperatures (above crossover) from a 
general master equation. Note that below the crossover 
the energy diffusive domain ceases to exist quickly with 
decreasing temperature. The main difference to previous 
treatments is that quantum effects like tunneling and re- 
flection are treated systematically and thus must also be 
taken into account in transition matrix elements, at least 
for energies near the barrier top. This way, we obtain a 
semiclassical generalization of the classical energy diffu- 
sion equation which will be used in this paper to analyze 
the impact of quantum fluctuations on the escape for a 
standard model of an energy barrier. In a subsequent 
publication the situation described above, namely, quan- 
tum effects in the switching between dynamical states 
will be treated. 

The paper is organized as follows: In Sec. |TT]we recall 
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the classical results and introduce the basic notations. 
In Sec. mil it is shown how to re-derive the classical en- 
ergy diffusion equation from a semiclassical limit to the 
quantum master equation, where tunneling is neglected. 
The latter process is taken into account systematically in 
Sec. IIV( where the main result, namely, a quantum ver- 
sion of the classical energy diffusion equation is derived. 
Quantum corrections to the escape rate for a standard 
potential model are obtained in Sec. |V] and discussed in 
Sec. Ell 



II. CLASSICAL REGIME 

In the standard description of dissipative systems the 
total Hamiltonian is of the form uQ, 



H = Hs+Hb+Hsb 
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where Hs contains the system part with a metastable 
potential 



90 



(2) 



with a well region around g = separated from a con- 
tinuum by an energy barrier of height AV = ^m^Q^Q. 
The heat bath Hb consists of a large number of harmonic 
oscillators linearly coupled via Hsb to the system. 

To setup the stage and for later purposes we first con- 
sider the well-known classical realm. There, after inte- 
grating out the bath degrees of freedom the associated 
equation of motion for the system coincides with a gen- 
eralized Langevin equation [26j 

dV /■* 

mq{t) + — +m ds^it-s)q{s)=C{t). (3) 
9q Jo 

Here the influence of the heat bath at temperature T — 
l/k^(3 is described by a damping kernel 7(i) which is 
related to the Gaussian noise force ^{t) with vanishing 
mean via the dissipation fluctuation relation 



t\ 



(4) 



In the sequel we will assume that the bath memory time 
is much shorter than other relevant time scales so that a 
Markovian kernel of the form ^{t) = 2'yS{t) applies. 

Thermally activated decay in the regime of weak dissi- 
pation requires that the energy lost by a classical particle 
during one period of oscillation in the metastable well is 
sufficiently smaller than the thermal energy. Specifically, 
as first shown by Kramers the condition reads 



-fI{E = AV) < 1//3 



(5) 



with the action I{E) along a periodic orbit being almost 
constant for energies close to the barrier top, which is 
the relevant energy range for the escape process. A con- 
venient procedure to derive the rate for this process con- 
verts the above Langevin dynamics in an equivalent equa- 
tion of motion for the probability distribution W{E, cj), t) 
to find orbits with energy E and phase (j) at time t in the 
well. Since for low friction energy is almost conserved, 
while the corresponding phase oscillates fast, the latter 
one can be adiabatically eliminated and one arrives at 
the energy-diffusion equation for the marginal probabil- 
ity P(i?,t) [l^ 



PiE,t) 



_d_ 

dE 



ME) 1 



l_d_ 

PdE 



uj{E) 
2ti 



P{E,t). (6) 



Here, A{E) is the energy relaxation coefficient 

A{E) = jI{E) . (7) 

The escape rate is determined by the stationary nonequi- 
librium distribution Pst{E) to ^ which is associated 
with a finite flux across the barrier and obeys the bound- 
ary conditions that Pgt = for E > AV and that Pgt (E) 
approaches a Boltzmann distribution in the well region. 
Accordingly, one obtains for high barriers /3AV ^ 1 the 
classical Kramers result 



dEP^tiE) 



uojHAVJl^ 
2tt 



-0AV 



(8) 



The above expression can also be derived by working in 
a multidimensional space including the system degree of 
freedom and the bath oscillators. One introduces normal 
mode coordinates in the parabolic range around the bar- 
rier top, where the total system is separable, and then 
studies the dynamics of the unstable normal mode in 
presence of the coupling to the stable ones due to the po- 
tential anharmonicity. This methodology allows to cap- 
ture the full turnover from the regime of weak dissipa- 
tion considered here to the regime of spatial diffusion for 
higher friction [13, ■ F'or weak dissipation one recovers 
to leading order ([8]), while an improved result including 
higher order corrections (but still away from the turnover 
region) reads Fd^imp — (Ab/wb) Tci- Here, finite recross- 
ing at the barrier top is described by the Grote-Hynes 
frequency [1^ 
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(9) 



which describes the effective barrier frequency of the un- 
stable normal mode, while Wb is the bare barrier fre- 
quency. For the model ^ one has uji, = ujo. 



3 



III. SEMICLASSICAL DERIVATION OF THE 
KRAMERS EQUATION 

In a first step and as test-bed for the semiclassical ap- 
proximation we derive here the classical Kramers equa- 
tion ([6|) from the Pauli-master equation. For this we fol- 
low the approach developed in [30| for potentials with 
stable bound states. Thus, we assume a number of 
N discrete energy levels in the well E < AV and at 
this point neglect quantum tunneling through the barrier 
completely (cf. fig. [1]). At the end, N will be taken to be 
very large so that one can work with a quasi-continuum 
of states. Accordingly, we have 



Hs\n)^E^\n), n = 1,2...N 
with the normalization 



(10) 



(11) 



Now, let Pnit) be the probability to find the system in 
discrete state n at time t. Then, one can derive from 
the exact path integral expression of the reduced density 
matrix for weak coupling to the heat bath and not too 
low temperatures Pauli's master equation, i.e., 



N 



(12) 



The transition rates from state m to n are given by 31 1 



dtTiB[{n\HsB{t)\m){m\HsB\n)p''B] 



(13) 

with Tr^ denoting the trace over the bath, the 

coupling in the interaction picture HsB{t) = 
e^(Hs+Hs)t/njj^^^-^KHs+Hs)t/n^ and 

g/3ffi3^^j.^jg/3ffi3j ^j^g gquilibrium bath density ma- 
trix. The above expression can be evaluated explicitly in 
case of the linear system-bath coupling and one arrives 
at the golden rule type of formula 



Wn,m = ^\{n\q\m)fD{En - E^) 

with the bath absorption/emission captured by 
D{E^ - Ek) = 2m-i{En - E^) n{En - Ek) , 
where 

n{E^ -Ek) ^ 



(15) 



(16) 



The probability Pn{t) is related to the probability 
P{E, t) to find the system in a state with energy between 
E and E + dE at time t via 



p(i;,t) = ^(5(£;-£;„)p„(t) 



(17) 



which allows to reformulate the above master equation 
(fT2ll as 



N 



P{E,t)^ J2 SiE-E,,)[Wn (t) - it)] 

m.n—O 
N N-n 

- ^ ^ [S{E - E„+i) - S{E - E,,)] Wn+l,nPn{t). (18) 



n=0 l=~r 



A. Semiclassical wave function 

To calculate the transition matrix elements in system 
space we now apply the semiclassical WKB approxima- 
tion [30|- Up to second order in h one obtains for the 
wave functions 

{n\q) EE {E,,\q) = 1 [(i?„|g)- + {E,M)+] (19) 



with 



N{E) 
VpiE^ 



(20) 



and the action So{E, q) = JJ^ piE, q')dq' of an orbit start- 
ing at qi and running in time t towards q with momentum 
p{E,q') at energy E |32]. The normalization is deter- 
mined from (jlip 



N{E) = y/2muj{E)/TT 



(21) 



with uj{E) being the frequency of the classical oscilla- 
tion at energy E. Note that forward and backward wave 
contributions in ([^0]) are related by 



[{E\q)+Y = {E\q)- 



(22) 



The energy En ~ E of the discrete spectrum is given by 
the quantization condition 



J^Idqp{E,q)^n+l 



(23) 



(14) containing the action over periodic orbits in the well. 



B. Semiclassical matrix elements 

With the semiclassical wave function at hand we can 
now calculate the transition matrix elements which enter 
the transition rates in (fT4|) . According to the restricted 
interference approximation [33] we only keep the diagonal 
contributions of forward/backward waves to obtain 



{E,n\q\En) 



dqq [{E„,\q)+ {En\q)^ + {E„,\q)- (E^lq)- 



n=0 



^-^jdq{E„Myq{En\q)-=Qtx 



1 ' 



(24) 
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where qi, (72 denote the left and the right turning points, 
respectively, of the periodic orbit with energy E. To 
calculate qIci™^ we exploit that due to [Hs, q] = hp/ (im) 
the matrix elements Q^"'™) = {Em\q\En) and p("^™) = 
{Em\p\En) are related via 



This way one finds together with (^0)1 



(25) 



-j(n,m) 
- scl 



h 1 



9p 



dq\ 



I uj{En)uo{E 

rn J 



2ip{E^)dq 



P{En)p{Ern) 
g^[So(B„,g)-So(i5,„,9)] (26) 



The exponential is further simplified by applying the 
semiclassical expansion of energy differences based on 
(ESI), i.e., 

hAi{E) En+i-En = niuj{En) +—P[Lo{Enf]'+o{K') , 

(27) 

where here and in the sequel the prime ' at en- 
ergy dependent functions denotes the derivative with 

respect to energy. Now, Qg[.{(-E„) = (3scl'"'^ ^ 
n = I can be systematically expanded in powers of h 
based on S'o(£'„,g) - So{Em,q) « -hluj{En)t{En,q) + 
^pLu{En){Lj{En)t{En,q)Y with t{E,q) being the time a 
segment of a periodic orbit with energy E needs to reach 
position q from its turning point. The expansion for the 
squared matrix element is thus found from (|26p to read 



(0^ 



(28) 



with 



(£;) = _L / dge-*'"(^'*(^-«) (29) 
2j7r{ J 

describing the classical transition amplitude. From this 
finding the semiclassical expansion of the transition rates 
(fH)) between states with energy En and energy Em such 
that En — Em — lhiij{En) is given by 



Wn,-^=Wl{En) = 

2m7 (;)2 IniLo-f 



hp 



(30) 



C. Semiclassical expansion of the master equation 

The above results are now used to determine the 
semiclassical expansion of the master equation ([TSl). 
which in turn means to treat a large number of quasi- 
stationary energy eigenstates in the well as a quasi- 
continuum. After some straightforward manipulations 



(see Appendix [X]) one arrives at 



P{E,t) 



^ ( dE 



k=l 
N 



d 



J2 [Wi{E)AnE) + W^i{E)Ali{E)] P{E,t) 

(31) 



1=1 



In case of a harmonic system in the second sum only 
terms with I = 1 contribute and together with the en- 
ergy independence of its frequency ujq the first sum gives 
rise to differential operators of the form exp{±fiLLjod/dE). 
For anharnionic systems we proceed in the spirit of the 
semiclassical expansion and take into account only terms 
up to fc = 2, i.e., 



PiE,t) 



(32) 



with 



N 

^^•(^)= E T^W,{E)Af{E). (33) 



Combining this result with pop provides the energy dif- 
fusion equation in leading order in H, namely. 



1 + 



l_d_ 

13 dE 



2-K 



P{E,t), 



where 



N 



A{E) = 2miu{E)j7T ^ Qi\\Ef 



(34) 



(35) 



In the limit of a quasi-continuum of states {N ^1) this 
sum can be recast with the help (|29p as 



A{E)^^ jdqp{E,q)^^I{E) (36) 

with the action of a periodic orbit with energy E 

I[E) = jdqp{E,q). (37) 

Accordingly, we recover from p4|) the classical Kramers 
equation dH), i.e. P{E,t) = Z:^'^ P{E) with 

=—^IiE)[l + -—^^. (38) 

The main result of this derivation is that up to second or- 
der in the ?i-expansion (|31|) quantum effects in the energy 
diffusion do not appear. Of course, the above procedure 
can be elaborated to systematically include higher order 
fi correction. We will see in the sequel though, that for 
the barrier crossing problem, finite reflection and trans- 
mission amplitudes provide much larger contributions, 
which are formally of order 1. 
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IV. ENERGY DIFFUSION FOR FINITE 
TRANSMISSION 

In the previous section the semiclassical wave func- 
tions Uni), (pOl) were taken as stationary eigenstates in 
the well potential for all energies below the barrier top. 
The underlying reasoning that finite barrier transmission 
does not play a role, however, fails as the temperature is 
lowered. Namely, for energies E < AV the finite proba- 
bility T{E) to timnel through the barrier and for energies 
E > AV the finite probability R{E) to be reflected from 
the barrier can only be neglected as long as near the 
barrier top the thermal energy scale ksT by far exceeds 
the quantum scale tuvh- When this does no longer ap- 
ply the consequences are two-fold: On the one hand in 
a range somewhat below AV discrete energy levels are 
strongly broadened due to tunneling so that one has a 
continuum of states in the well strongly coupled to the 
continuum on the right side of the barrier (see figd]); on 
the other hand, continuum eigenstates with E > AV gain 
a longer life-time in the well region due a finite reflection 
and may thus influence the escape process as well. Note 
that corresponding quantum effects are substantial since 
R{AV) = T{AV) = 1/2, and thus dominate against 
higher order ?i-corrections (associated with higher than 
second order derivatives in E) in the expansion (pij) [sot . 
However, it is well-known that in the limit H —^ the 
functions R and T are non-analytical in h, [see (1551) ] so 
that an expansion in powers of h is not feasible. In a 
semiclassical approach we thus proceed to work in a con- 
tinuum representation right from the beginning and treat 
low lying states and higher lying states on equal footing. 
This approach is justified as long as temperature is not 
too low and tunneling prevails near the top of the barrier. 
Hence, instead of (fT7|) we now start with 



Then, R{E) pE{t)dE is the probability to find the system 
at time t in a state with an energy between E and E + dE 
inside the well and T[E)pE{t) with T{E) = 1 - R{E) 
the probability for tunneling out of the well. In case 
of finite barrier transmission apart from transitions for 
states inside well induced by the heat bath, there is also 
a loss of population due to tunneling through the barrier. 
Hence, the variation in the probability psit) to find the 
system in a state with energy E in the well is equal to the 
probability that the system arrives in the well from an 
other energy E' minus the probability that the system 
leaves the state with energy E to another state in the 
well region and minus the probability of tunneling out of 
the well, i.e., 



PEit) 



dE' 
We'e 



We, 



R{E')pE'{t) 



n(E' 

R{E)pE{t) 

n{E) 



(40) 



nEf-^PE{ti 



where n{E) is the density of states. Now, in generaliza- 
tion of we have 



P{E,t) 



dE'dE" 



We" 



R{E')pE'{t) 



n{E') 



-We' 



R{E")pE"{t) 
n{E") 



S{E ~ E") 



-T{E)^^PE{t) 



(41) 



with We,e' being the transition rate from state with en- 
ergy E to state with energy E' , namely. 



E,E' 



1 f°° 

-2 / dtTrB[{E\HsB{t)\E'){E'\HsB\E)p%'] 



-\Qq^{E',E)\^D{E^E'), 



(42) 



V 




FIG. 1; Typical metastable well potential discussed in the 
text with quasi-stationary states deep in the well and a contin- 
uum of states near the barrier top, around which a parabolic 
barrier approximation applies. 



PiE,t) 



PE'{t)6{E~ E')dE'. 



(39) 



where Qq^{E' , E) = {E'\q\E). Here, the bath spectral 
function D{E) is defined according to ^T5\i . The semi- 
classical wave functions entering the transition matrix 
elements for the system degree of freedom, however, must 
now include the finite transmission/reflection probability. 
In the energy range close to the barrier top, where tun- 
neling dominates in the temperature range considered, 
the WKB approximation is not applicable because the 
classical turning points to the left and to the right of the 
barrier are not sufficiently separated. In this situation, 
one exploits that the Schrodinger equation for a parabolic 
barrier with barrier frequency tOh can be solved exactly. 
The proper eigenfunctions are then matched asymptoti- 
cally (sufficiently away from the barrier top) onto WKB 
wave functions to determine phases and amplitudes of 
the latter ones. This way, one obtains 

{E\q) ^NiE) [{E\q)- + riE){E\q)+] (43) 

with the matrix elements (j20p and the complex valued 
reflection amplitude r{E) of a parabolic barrier [s^l re- 
lated to the reflection probabihty R{E) = \r{E)\'^. The 
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normalization follows from {E\E') = 5{E — E') as 



NiE) 



/^V l + RiE) 



(44) 



In case of vanishing transmission, E « AV^ and i? ^ 1, 
one recovers the previous result N{E) so that 

it is possible to use p3|) for all energies provided the 
length scale where a parabolic approximation for the bar- 
rier applies is much larger than the quantum mechani- 
cal length scale i/sTmwb. This in turn requires a high 
barrier and thus follows directly from the condition for 
metastability. 



Semiclassical matrix elements 



In a similar way as in Sec. IIIIBI we now calculate 
in the semiclassical limit the transition matrix elements 
QqmiE' , E) for finite transmission through the barrier. 
One has 

Qq^{E',E) ^ ^(E')^(E) r dq[(E'\q) + q{E\q)- 



+r{E)r{Er{E'\q)-q{E\q)+] 
N{E')N{E) [Q:,i + r{E)r{ErQsci] 



where in generalization of ([24|) 

f92 



1 

Qscl{E',E)^- dq{E'\q)-q{E\q)- 



(45) 



(46) 



Note that in contrast to Q, 



(n,m) 
scl ' 



which is real, the tran- 



sition element Qqm is complex for finite tunneling ampli- 
tudes \r\ > 0. Hence, Qqm cannot be written as a loop 
integral as in To proceed, one introduces 



(47) 



E' -E = 6. 
Now, exploiting again ([25]) one has 



(S) 



where 



and 



huj{E) 
2m 5 



dq e 



K 



-it(E,q)S/h 



h dp{q) 



-it{E,q)S/h 



P'il) 
Pil) 



(48) 



(49) 



(50) 



ip{qY S dq 



92 - gie-^"*/'"'"^^) 



This latter contribution is a boundary term which van- 
ishes for the closed integral in The ?i-expansion of 
the squared matrix element thus reads 



\Q 



qm| 



qin 



(51) 



with coefficients Aqm, Bqm, and Cqm specified in Ap- 
pendix |B] 



B. Expansion of the master equation 

We now proceed along the same way as done for a 
discrete spectrum in Sec. pil Cp with the substitutions 
hAi S and J^i ^ ^ j Thus, we obtain from 
the expansion 



P{E,t) 



k=l 



_d_ 

'dE 



1 
fc! 



R{E)P{E,t) 
n{E) 



d5Ws{E) {-5)'^ 



T{E)^P{E,t) (52) 



with Ws{E) = We,e' for E' - E = 6. The leading order 
terms in the sum above with fc = 1 and k — 2 are kept to 
get the energy diffusion equation for finite transmission 
in the semiclassical limit, i.e. P{E,T) = d"^^'^ P{E,t), 
where 



C 



(scl) 



d_ 
'dE 



d_ 

'dE 



R{E)-T{E) 



u}{E) 
2ti 



Here, the moments of the energy fluctuations read 



(53) 



(54) 



1 /i" 



Apparently, when comparing this expression with (j33[ 
one observes that drift and diffusion coefficients Dk cor- 
respond, as expected, to {S'') — {{E' — E)'') for fc = 1, 2 
respectively. Now, combining the expansion for the den- 
sity of states 



n(E) 



1 



huj{E) 



0{n) 



(55) 



with Eqs. (P5)) - (|5ip one finds upon evaluating the mo- 
ments (j54p the first main result, namely, the semiclassi- 
cal expression (for details see Appendix [C| of the evo- 
lution operator in the energy diffusive regime P{E,t) = 
4^V(£;,i),i.e., 



,(scl) 



_d_ 

'dE 



C{E)^I{E) 1 



ld_ 

pdE 



uj{E)R{E) 



2-K 



T{E) 



2tt 



with 



C{E) 



, 1 + R{Ef 
\l + R{E)Y 



(56) 



(57) 
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Some remarks are in order here: First, for vanishing 
transmission {R — 1, T = 0) one recovers from the above 
expression the classical diffusion operator ([38|) . This in 
turn proves what we have said above, namely, that on 
this level of semiclassical expansion an energy diffusion 
operator can be derived starting either from a discrete 
or a continuous spectrum in the well. The difference is 
though, that the latter procedure conveniently accounts 
for barrier tunneling near the barrier top with the prop- 
erty that C^^'^^^ C^^'^ in the range E <^ AV, where 
R{E) 1. Second, higher order h corrections in the ex- 
pansion (|52[) can now be calculated accordingly, where, 
however, in case of a barrier crossing problem correspond- 
ing contributions are much smaller than those originating 
from quantum transmission and reflection. Third, in var- 
ious previous works an operator similar to the above one 
but with the factor C{E) = 1 has been used to study 
quantum effects in thermal activation in the low damp- 
ing regime (20l [2ll . [2^ . [ssj . The important difference is 
though that here we derived this operator including its 
transition matrix elements within a consistent semiclas- 
sical expansion of the underlying wave functions, while 
in the former cases in (|4ip the corresponding classical 
expressions were adopted. Hence, the above equation is 
the systematic generalization of the classical Kramers' 
equation in energy space to lower temperatures. 



V. QUANTUM ESCAPE RATE 

The escape rate in the quantum regime can now be 
evaluated similar as in the classical case by searching for 
the quasi-stationary energy distribution Pgt {E) . For this 
purpose, we start by specifying the known transmission 
and reflection probabilities in a uniform semiclassical ap- 
proximation, i.e.. 



T{E)^\t{E)\' 



1 



l+CXp[-S'e(^)/?i] 



(58) 



where Sc{E) denotes the Euclidian action of a periodic 
orbit with energy E oscillating in the inverted barrier 
potential —V{q). In the energy range near the bar- 
rier top, where tunneling dominates in the temperatures 
range considered here, the action reduces to its vale for 
a parabolic barrier, namely, 5'e,pb(-£') = 2Tr{E — AV)/ujh- 
The above expressions thus smoothly connect the energy 
range near top with the low energy range, where T{E) 
drops exponentially. Hence, deep inside the well the only 
quantum effects are zero-point fluctuations. The reflec- 
tion dependent factor C{E) is always larger than 1 which 
means that leaking out of the wave function by tun- 
neling appears effectively as an increase in energy loss 
C{E)jI{E) > ■yI{E) during one cycle of a classical orbit 
in the well. 

Now, in order to explicitly find the steady-state dis- 
tribution Pst{E) it is convenient to work with a di- 



mensionless energy measured from the barrier top, e = 
{E—AV)P, and to introduce a dimensionless inverse tem- 
perature 



(59) 



Then, one writes /(e) = [tj(e)/27r] i?(e)Pst(e) so that from 
Pst — one arrives with ((56l) at 



(60) 



The boundary conditions here are such that /(e) — > 
for e ^ 1 and that inside the well (e some 9 below AV^) 
/(e) — !■ /^(e) with the equilibrium distribution 



/Me) 



1 



2TThZf) 



(61) 



Here, the quantum partition function in the harmonic 
well reads 



Zn = 



1 



n 



n=l ^" + ^0 + '^nl 



(62) 



with Matsubara frequencies Vn = 27m/?i/3. For van- 
ishing friction this reduces to the known result Zqq — 
l/[2sinh(wo?i/3/2)]. With 



/(e) = Me) g(e) 



([50)1 becomes 



C{e)g{e)"+g{ey[C{ey-C{e)] 



lT(e) 



pR{e) 



(63) 



5(e) , (64) 



where we also exploited that the action varies smoothly 
in the relevant energy range around e = (i.e. in an 
interval of some k^T^ fiujiy below the top) so that with 
/(e) ~ /(e = 0) one defines p = PjIiO). 

Analytical progress can now only be made by approxi- 
mating the coefficient C(e). For this purpose we assume 
self-consistently that g{e) becomes exponentially small in 
the semiclassical sense for energies e of order 6 above the 
barrier top. Then, fM]) needs only to be considered in a 
relevant range e < 6- k with some constant k > 1 of order 
1. A very accurate approximation in this energy range is 
given by 



c{e < e ■ k) K 1 + 



(65) 



for K < 1.5. Since the final rate does not depend on the 
precise value of n, we choose k = ln(4) which ensures 
that C(e)' is very well approximated for all energies up 
to its maximum located at e = ln(4) 6. The solution of 
(j64p is then found in terms of hypergeometric functions 
as 



Be' 2F1 



a, 1 



e^ 
"9" 



--a,- + -+a,l 



(66) 



8 



with the coefBcients 



B = 



- a, ■ 



0,1 + 6*, 



and 



- 6')2 + 366*2 
4^ 



(67) 



(68) 



A straightforward analysis now verifies our assumption 
that g{e) ^ Q for e of order 6 above the barrier top. 

The semiclassical approximation to the quantum es- 
cape rate 



r = 



' dE T{E)'^P,,{E) 



27r 



(69) 



is obtained by exploiting Pst — with ([SS]) and (ISS)) 



as 



(70) 



Here, the lower limit lies in an energy interval some 9 
below the barrier top where the nonequilibrium distri- 
bution matches onto the Boltzmann distribution in the 
well. In this range C(e) ~ 1, while still /(e) w I{e = 0) = 
I{E = AV). To leading order in the friction strength we 
thus find the second main result, i.e., 



_ sinh(u;on/3/2) 



(71) 



with the classical result specified in ([8]). The first factor 
captures quantum effects (zero-point fluctuations) in the 
well distribution, while the second one describes the im- 
pact of finite barrier transmission close to the top. More 
details will be discussed in the next Section. 



VI. DISCUSSION 

Of particular interest are the leading quantum correc- 
tions to the classical rate expression when either of the 
parameters p or 6 becomes small. For somewhat larger 
friction (but still sufficiently away from the turnover 
range) higher order corrections can also be accounted for 
along the lines described in f25| . 



A. Quantum effects 

To further analyze quantum effects in (|7ip we expand 
the prefactor y = Tsd/^d = y{&,p) in the inverse tem- 
perature 9 as well as in the friction strength p. From 
one observes, however, that the limits p and 9—^0 
are not interchangeable: In the classical range and for 



£5 




FIG. 2: Escape rate normalized to the classical rate (|8} versus 
the dimensionless inverse temperature 6 — 'hLU\^j3. Thin lines 
refer to the result from [U] where C{E) = 1 [cf. [56])], thick 
lines to the new expression (|7ip for 'y/uoo =0.0015 (solid), 
0.005 (dashed), 0.01 (dotted). 

energies below the barrier top the right hand side van- 
ishes according to T{e)/R{e) = exp(e/6') 0, while for 
any finite 9 it diverges for p ^ 0. This behavior refiects 
the fact that quantum mechanically a steady state with 
a finite flux out of the well exists even in absence of dis- 
sipation. Hence, we consider flrst for flxed p quantum 
corrections to the classical high temperature limit 9^0 
and gain 



y^l-9bi+9^{b2 + TT^/G) + 0{9^/p) 



(72) 



Here, the coefficient bi ~ 1.04 originates merely from 
the expansion of B and describes the impact of finite 
barrier transmission, while 62 ~ 0.20 is determined also 
by the well partition function. An inverse friction depen- 
dence appears in third order terms, which in turn means 
that the above expansion is only valid provided p remains 
finite and that its range of validity shrinks with decreas- 
ing p. Interestingly, there exists a temperature range 
with an upper bound given by 61/(62 + ""^/B) ~ 0.56, 
where quantum effects reduce the escape rate below its 
classical value and only for lower temperatures, i.e. larger 
6*, does the quantum result exceed the classical one (see 
fig-El)- This finding has already been predicted in [2lll2^ 
with a quadratic temperature dependence in leading or- 
der though, which relates directly to the fact that clas- 
sical transition matrix elements have been used in these 
previous works. Apparently, the impact of quantum fluc- 
tuations in (|7ip is enhanced and the suppression of the 
escape due to a finite reflection at the barrier top is 
diminished. It is only for lower temperatures that the 
quadratic dependence on the inverse temperature dom- 
inates and always leads to a rate increase compared to 
the classical result. 

The leading quantum correction for p — > at flxed 
temperature can already be read off ([M)) . Namely, the 
energy range below the barrier top where the flux solution 



9 




-4 -3 -2 

Logio(7/'^o) 



must be used in the well together with a factor describ- 
ing the influence of the stable normal modes around the 
barrier top. This way, one gets 

1 °^ 2 1 
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FIG. 3: Escape rate normalized to the classical TST rate 
(Ftst ~ ijJo/2-Ke~^^^) as a function of the dimensionless 
friction strength for various values of the dimensionless inverse 
temperature 9 — hLJhl3/2n {9 = denotes the classical rate) 

matches onto the thermal equilibrium grows for p ac- 
cording to e ~ 0lii{p). To leading order one then obtains 
in this domain d^g oc — exp(e)/p^, where the proportion- 
ality constant is determined by the boundary conditions 
and must be obtained from the full solution The 
leading order quantum correction to the rate prefactor 
scales thus like y cx , which is non-analytical in p 
[2^ . The detailed calculation using large parameter ex- 
pansions of the hypergeometric functions [36] gives 



3^^ 



sinh(wo?i/3/2) F(l - 9) 



-0{p 



-e+i/2\ 



(73) 



This behavior is illustrated in fig. O where the rate be- 
comes extremely sensitive to 9 for weaker friction and 
lower temperatures. It is interesting to note that for 
p and sufficiently high temperatures our result ([73]) 
coincides with the result obtained previously in (2l| since 
in this limit the matching happens for energies below the 
barrier top where C = 1 in (j64p . 



B. Higher order friction corrections 

As already mentioned at the end of Sec. [ll] the energy 
diffusive limit can also be described by a normal mode 
analysis in full configuration space. In [22| this approach 
has been applied to include quantum effects due to a 
finite barrier transmission, while transition matrix ele- 
ments have been treated purely classical. Based on the 
above analysis the generalization of the latter ones to 
the semiclassical domain is straightforward. Of course, 
to leading order in the dissipation strength one recov- 
ers ([7T|) . Higher order corrections (still above turnover 
though) appear on the one hand through cjb — ^ Ab [see 
^] in the factor B, thus capturing the impact of friction 
in the tunneling process. On the other hand, the full par- 
tition function (|62p of the damped harmonic oscillator 



APPENDIX A: MOMENTS OF ENERGY 
FLUCTUATIONS 

In this Appendix, we specify how to write the master 
equation (jl8|) in a way, which is amenable to a system- 
atic semiclassical expansion. Therefore we can use the 
?i-expansion of the delta function 

5{E - En - hAi{En))f(E,,) = [5{E - S„) + 



E 

fe=i 



d_ 

dE 



{-hAijEr,))' 

k\ 



■S{E-Er,)]fiE„ 



and rearrange the sum in (|18p in the following way. 



N N-n 



N 



n— l = -n 
1^0 



1 = 1 



N 



^(/(n,0 + /(n,-0) 



.71=0 



l-l N 

E/("'-o- E /("'O 

n=N-l+l 



n=0 



In (UHl), there is f{72,l) = Wi{En)Pn and 

W-iiEn) - {En-i\q\Er,) =0 ff n < / - 1 
Wi{En) =^ {En+i\q\En) = if n > N - I + I 

Therefore the last two sums in (|Aip vanish. Eqs. (|Aip 
and ljAl[) may be inserted in (fTS]) to give the final result 
given in (1311). 



APPENDIX B: COEFFICIENTS FOR 
TRANSITION MATRIX ELEMENTS 



Here we collect the coefficients appearing in the h- 
2 



expansion of [Q^ml^ in ([?T|) . One has 



iqm = |Ocl|'(i?' + 1) + RiQll + Qll ) 



and 



Bq„, = 2iV37V'[|Qei|'(i?' + 1) + R{Ql + Qli)] 
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as well as 

+R\Qci\^R' + r*r'Q*l + rr*' Q\ 

+ {R' + l){QdK* + Ql,K) + 2R{Q,iK + Ql,K*) . 



APPENDIX C: EXPANSION OF THE 
DIFFUSION OPERATOR 

In this Appendix we will expand the moments of the 
energy fluctuations {5) and (5^) in order to be able to 
write a semiclassical expression for /I^'^'' . 

In a first step we write {5) with US]) , (|42]) , ^ and ^ 

as 

With the expansion of the Bose occupation function 

1 



n{5) 



_ 1 1 

e'^* - 1 ^ J5^2 



with ([5T|) and the results of IB] we obtain the ?i— expansion 



d55' 



(36 



+ 



{B + N^C) 



The integration leads to 



(-5) 



27r 



+ l)iv37VV/ + l^-^7V4(cj2/)' 



, (CI) 



where J denotes the action (157]) . 
In the same way one arrives at 



ijj{E)m'y 



..N-^A 



27r/3 



+ 1) 



(C2) 



Now we are able to combine (|Cip and (|C2[) with ([55)1 to 
get the evolution operator in the energy diffusive regime 
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